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, The two-fermion system covers a huge number of applications in atomic (e + c ), nuclear (NN, NN) and subnuclear 
(qq) physics. The interest in using a relativistic description for such systems appeared in the early days of quantum 
mechanics 1,;2] and has constantly been pursued since by manyauthors. This interest has recently found a new elan 
due to the measurements performed at Jefferson Laboratory 0,0, H @ where simple nuclear systems have been - and 
are being - probed at momentum transfers much larger than their constituent masses. This experimental activity 
. motivated a consequent number of works on relativistic dynamics. Extensive reviews on the past and recent deuteron 
Q\ ' results can be found in 0, @ • 

Most of the approaches developed for describing relativistic two-bod y sy stems are based on Bethe-Salpeter (BS) 
equation [![Ii[limHQ[liIor three-dimensional reductions of it Hi E3 E H H E3 . 

An alternative approach is provided by the Light-Front Dynamics (LFD). In its standard version, following Dirac's 
classification of relativistic theories [22} , the state vector is defined on the a = z + t surface. Wave functions - defined 
as the Fock components of the state vector - are the formal objects of this theory and are directly comparable to 
their non relativistic counterparts. LFD has been developed and used b y many authors jl^, 03 IH Eo, 113, |2^| EH 
IMISSHIMIMIIllMlMIIilMH represents a promising approach 

to non perturbative Hamiltonian Quantum Field Theory, specially when dealing with composite relativistic systems. 
The interested reader can be aware of the last advances and more complete references set in the proceedings of the 
r-j last conferences devoted to the subject [52ll53| . 

The explicitly covariant version of Light-Front Dynamics (ECLFD) was initiated by one of the authors in a series of 
. £h ' papers |54l l55l l56| . The state vector is there defined on a space-time hyperplane whose equation is given by ui ■ x = a, 
where w is a four-vector determining the orientation of the light-front plane and satisfies lo 2 — 0. This choice is not 
only a mathematical delicatesse but a way to carry everywhere in the theory the w-dependence in an explicit way. It 
has several advantages, all related to the fact that a; is a four vector with well defined transformation properties. This 
approach provides explicitly covariant expressions for the on shell amplitudes, a property which is often hidden in the 
standard formulation, recovered by fixing the value ui — (1, 0, 0, —1). This value is however associated to a particular 
reference frame and it is not valid in any other one. The formalism and some of its first applications to few-body 
systems has been reviewed in |57^ . 

Approximate light-front solutions for the NN system |58l l59j were found in a perturbative way over the Bonn model 
wave functions [6(| and successfully applied to calculate the deuteron electromagnetic form factors [6l| measured at 
Jefferson Lab. Latter applications to heavier nuclei [62l Ir33l ] have shown the pertinence of this approach in describing 
high momentum components of the NN correlation functions. 

These successes stimulated a series of works aiming at developing some formal problems of the theory and to obtain 
exact solutions in the ladder approximation for systems of increasing complexity. Results concerning bound states of 
two scalar particles can be found in |6J,[65, 66, 67]. We present in this paper the formalism and numerical solutions 
describing bound two fermion systems interacting via the usual - scalar, pseudoscalar, vector and pseudovector - 
one-boson exchange (OBE) kernels. Results are limited to J = and J = 1 states. Our main interest in this work is 
to study the solutions of the LFD equations as they are provided by the OBE ladder sum with special interest in their 
stability, their comparison to the non relativistic limits and the construction of non-zero angular momentum states. 
For this purpose, we have studied each coupling separately and the only physical system considered is positronium. 
The first conclusions concerning the Yukawa model have been published in (6^, kXI. Fol Flj and a more detailed 
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derivation of equations and kernels can be found in |72| . This series of works is also being extended to the two-body 
scattering solutions and to three-particle systems. The case of three-bosons interacting via zero range forces was 
considered in |23|. In references the ensemble of these results is briefly reviewed. 

It worth mentioning preceding works on two-fermion system using the LFD approach. In |46|. the relativistic 
bound-state problem in the light-front Yukawa model was considered. In |37l l38| , positronium and heavy quarkonia 
calculations in discretized light cone quantization were carried out. The formalism was used in |27j ] to build one boson 
exchange kernels and to calculate nucleon-nucleon phase shifts as well as deuteron properties. Recent application to 
meson spectra can be found in |40t |4l| . LFD was also applied in |43J, |45( to describe the NN system and nuclear 
matter equation of state. 

The paper is organized as follow. In section [n] we establish the structure and main properties of the explicitly 
covariant Light-Front wave functions, the two-body equation and the OBE kernels. In section ITTT1 the problem of 
angular momentum J is discussed and states with J = 0,1 are constructed. In section llVl we derive the coupled 
equations for the wave function components of states with angular momentum J = 0. The corresponding equations 
for J = 1 states are derived in section The non-relativistic limit and perturbative calculations are discussed in 
section fvTI In sections IVIII IVlfTI and Hxl we present the results of numerical calculations. In order to disentangle 
their different behaviors each coupling is separately analyzed. Section Q contains a summary of the results and the 
concluding remarks. 



II. WAVE FUNCTION, EQUATION AND KERNELS 




FIG. 1: Graphical representation of the Light Front two-body wave function. Dash line corresponds to the spurion (see text). 



Wave functions we deal with are Fock components of the state vector defined on the light-front plane lo-x = 0. For 
a two-fermion system - shown graphically in Fig. it reads: 

were <7j are the constituent angular momenta. The general form of the wave function is obtained by constructing all 
possible spin structures compatible with the quantum numbers of the state. The four-vector u> enters in the wave 
function on the same ground than the particles four-momenta, giving rise to a number of structures larger than in 
non relativistic dynamics. Each of them is mastered by a scalar function, denoted fi all through the paper, which 
can be interpreted as a wave function component on the spin space. The number N of such independent amplitudes 
simply follows from the dimension of the spin matrix forming the two-fermion wave function with total momentum 
J, i.e. N = \(2J + 1)(2ct 1 + 1)(2ct2 + I) with a factor i to take into account the parity conservation. In the case 
ci = o"2 = h, it gives N=2 amplitudes for J=0 states and iV=6 for J=l. These wave function components will be 
specified in the coming sections. 

Since the Fock-space component is, by construction, the coefficient of the state vector decomposition in the cre- 
ation operators basis: aj^ (fo)'^ (/ci)|0), the independent variables are the three-dimensional vectors (k\,k2) and the 
particles energies are expressed through them. Consequently all four-momenta are on corresponding mass shells: 
kf = k 2 = to 2 , p 2 = M 2 , (lot) 2 = and satisfy the conservation law: 



k\ + ki = p + LOT. 



(2) 
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This equation generalizes the (_L, +)-components conservation in the standard approach; the minus components 
are not constrained. In the light-front coordinates with lo — (1,0,0,-1), the only non-zero component of lo is 
c<j_ = loq — w z = 2. The four- vector lot just incorporates the non- vanishing difference 2r = + k%- — P-. In this 
sense the ECLFD wave function is off energy shell. Since the four-momentum lot enters in the wave function on equal 
ground with the particle momenta, we associate it for convenience with a fictitious particle - called spurion - showed 
in Fig. Qby a dash line. We would like to emphasize however that the Fock space basis does not contain for all 
that any additional and unphysical degree of freedom. By spurion, we mean only the difference - proportional to lo 
- between non-conserved particle four-momenta in the off-energy-shell states. 

It is convenient to introduce other kinematical variables, constructed from the initial four-momenta as follows: 



k = L- x {V)k x = jfei - 
n \L-\V)Q\ 



V 



'V 2 + T 



(3) 



where V = p + lot, and L 1 (P) results from the Lorentz boost into the reference system where V = 0. In these 
variables the wave function JIJ is represented as: 

£ CT2CT1 =<P CT2CT1 (M). (4) 

Under rotations and Lorentz transformations of four-momenta ki, k2,p, lot, variables (k,n) are only rotated, so the 
three-dimensional parametrization Q is also explicitly covariant. In practice, instead of the formal transformations 
(|3*|) . it is enough to consider the wave function and the equation in the cm. system where V = k\ + ki = and 
set k\ = k, = —k, lo = fi\C3\. Because of covariance, the result is the same as after transformation J2J. Since 
lo determines only the orientation of the light-front plane, the modulus \C3\ disappears from the wave functions and 
amplitudes. Note that in the cm. system, the momentum p is not zero: p — —lot. 

The light-front graph techniques is a covariant generalization of the old fashioned perturbation theory. The latter 
was developed by Kadyshevsky [?7j and adapted to the explicitly covariant version in p3 • 

The equation for the wave function is shown graphically in FigEl It is the projection on the two-body sector of the 




FIG. 2: Equation for the two-body wave function. 



general mass equation P 2 (f> — M 2 4>. Its analytical form is obtained by applying the rules of the graph techniques to 
the diagrams in Fig. In variables @ this equations reads: 



4(fc 2 +m 2 ) -M 2 ] £ CT2CT1 (M) = -£3 / J2 lCthk,k',n,M 2 )<P a , a ,(k', 



d 3 k' 
) 

Efc' 



(5) 



where Kalal{k, k' ,n, M 2 ) is the interaction kernel. We detail in what follows the LFD one-boson exchange kernels 
corresponding to the interaction Lagrangians: 

(t) Scalar (S): 



(6) 
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(ii) Pseudoscalar (PS): 



(Hi) Pseudovector (PV): 



(iv) Vector (V): 



with 



C int = i g ps ^ (ps) (7) 



C mt ^7 M 7s 1> ^ (ps) (8) 



C mt = ^[g v r^ ] + ^ v {0^ ~ (9) 



fir I 

a = 2 w 7 - 7 7 )• 



W\ 2 k 1 ! k 
— > ^ > — /c 



7 co x 



COT <-> / CO^ - s 7 



k'„ 1 K 



2 K 2 ^ rv 2 

FIG. 3: One boson exchange kernel. 

The LFD ladder kernels have two contributions corresponding to the two time-ordered diagrams (in the light-front 
time) shown in Fig. |3| For S, PS and PV couplings they have the structure: 

X^^(*i»*a.^;*i»*2.^) = -T^2 W 2 {k 2 )0 2 u^{k' 2 )\ 1^(^)0^ (k[)} (10) 

4m L J L J 

6 (u-(k[ - h)) 6 (u-(ki - k[)) 



ji 2 - (k[- h) 2 + 2rV(fc 1 - ki) /i 2 - (*i ~ K? + 2ro;-(fci - *i) J ' 
For scalar exchange 

Oi = 2 =g s , 



for pseudoscalar 



and for pseudovector 



Ot 



O 



2 



Oi = 2 = H59p S 



(i + ^ns/^, if w(fei-*i)<o 

( 1 + f^) l 75/ pu , if co-fa -fcj)>0 
(i-fe)^/^, if w(*i-Ai)<0 
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with: 



4e 2 - M 2 , 4ei - M 2 

T= , T — 



2uj-p 



2uj-p 



For values r, t' ^ the kernels are off energy shell. In this case the pseudo-scalar and pseudovector kernels differ 
from each other but coincide on energy shell (r = r' = 0). 

We use the notation w — w M 7 A1 . Writing the propagators in the center of mass variables, (|1U|) gets the simpler form: 



l <72fl 



1 



1 



4to 2 Q 2 + 



[u„ 2 (k 2 )0 2 u a > 2 {k' 2 )\ [n CTl (fci)OiU^(fci)] , 



with 



{k-k'Y - (k-n)(k'-n) 



k-n k'-n 



(11) 



(12) 



The kernel for the vector coupling is given by a contraction of terms similar to I jllj l with the tensor structures L a p. 
It reads 



1 



1 



4m 2 \i 2 + Q 



:L al3 [u(h)Ofu(k[)} u{k 2 )O 2 u{k' 2 ) 



with 

L a p — 
and vertex operators 

Of 



-9a/3 + ^(fa - K - uT) a (k 2 - k 2 - ljt')p, if u)-(k x - k[) > 
-g a p + -p-iK - h - UT') a (k 2 - k' 2 - ujt)p, if w-(fei-fc' 1 )<0 



S»7 a + 4^' a (-«)(fci - *1 - urdaf, if - fci) > 

gvl a + 4^v ala {i)(k' 1 -k 1 -uT , ) a ,, if «■(*!- A^)<0 



(13) 



(14) 



(15) 



a, 



i-lot')p>, if uj-{k\ - k[) > 
+ JX'^HX^ -k' 2 - wr)p, if w(fc x - < 



(16) 



Hereafter we will not take into account the tensor coupling, that is we put f t = and Of = 2 = g v J a - In this 
case, vector kernel (|13|) simplifies into: 



I(- Cr 2 CT l 



9v 



Am 2 \x 2 + Q 2 



[tZ(fci) 7 Q u(*;i)] [«(*2)7a«(^ 2 )] 



[w(fci)u>u(fcj )] [u{k 2 )uju{k 2 )] 



(17) 



In the fi = case, e.g. one-photon or one-gluon exchange kernels, the L a p expressions depend on the gauge. Using 
the Feynman gauge, one has L a p — —g a p, i.e. the /i-dependent terms on (|14|l and (|17JI are simply dropped out. 

It will be often necessary to regularize the LFD kernels by means of vertex form factors. Unless explicit mention of 
the contrary we will take the form used in the Bonn model |60j| . i.e. 



F(Q 2 



A 2 - [i 2 
A 2 + Q 2 



(18) 



where A and n arc parameters whose values depend on the coupling. Form factors appear in the kernels multiplying 
each of the vertex operators Oi. In the non relativistic limit, Q 2 « (k — k') 2 and F is local in configuration space. 
This locality is however broken from the very beginning in LFD due to the n-dependent terms on Q 2 . 
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III. ANGULAR MOMENTUM 



In LFD the construction of states with definite angular momentum is a delicate problem. Working in the explicitly 
covariant version, we have developed a method to overcome this difficulty. It will be explained in this section. In 
contrast to the equal-time approach, the LFD generators J pl = J° 7 + J™* of four-dimensional rotations are not 
kinematical, but contain interaction in J™*. The interaction also enters in the angular momentum operator, i.e. the 
Pauli-Lubansky vector: 



Like the action of the Hamiltonian on the Schrodinger wave function is expressed through the time derivative 

H mt </> = id t 4> , 

the action of J*™' on the LFD state vector is expressed through derivatives with respect to the four- vector u> : 

j™*0H = v( w )0H, (20) 

where: 

d d 



Equation (|2U|I is called angular condition and can also be written in the form: 

</>(w) = Wf, </>(u) (22) 

with 



W M = -e^P'M^ (23) 



and 

Wfj, is a kinematical Pauli-Lubansky vector. As far as the angular condition is satisfied, the dynamical Pauli-Lubansky 
vector can be replaced by the kinematical one W^. The great benefit of doing so is that the problem of constructing 
angular momentum states with operator (|23[) becomes purely kinematical. In practice one rather prefers to start 
constructing states with definite angular momentum using Wu, and then take into account the restriction imposed 
by the angular condition 12U|) . 

It is worth noticing that without this condition, there is an ambiguity in defining the state vector with given angular 
momentum. This can be seen by introducing the operator: 

*-m'- 

It commutes with P M and W M and - taking A 2 instead of A - with the parity operator. The state vector is then 
characterized not only by its mass M 2 , momentum p, angular momentum J - defined by means of i|23fl - and parity 
7r but also by a, square root of the A 2 eigenvalue: 

A 2 4>^ = a 2 4> {a) . (25) 

For a total angular momentum J there are J + 1 eigenstates cf>( a \ In principle one could imagine any of these 
eigenstates to be an acceptable solution. It turns out however that, except for J = 0, none of these eigenstates can 
satisfy the angular condition l|22l) . Indeed if 4>(oj) is an eigenstate of A 2 , the right hand side of l)22|l - W^^ui) - is 
still an eigenstate of A 2 whereas this is not possible in its left hand side - S^(j>(uj) - due to the non zero commutator 
[Sp, A 2 ] 7^ 0. What is then the state vector? 



7 



A solution of the angular condition - the only remaining equation to be fulfilled - is therefore provided by a linear 
combination of different eigenstates : 

J 

a=0 

The coefficients c a can in principle be determined by inserting l|26[) into l|20|) or (|22ll . 

We would like to emphasize this result, which is, to our opinion, an important issue of Light-Front Dynamics. It 
tells us that the state vector is necessarily a superposition of different A 2 eigenstates. This conclusion does not depend 
on the approximation resulting from any eventual Fock-space truncation. 

In an exact solution of the problem, i.e. with the generators satisfying the Poincare algebra, the eigenstates 4>^ are 
degenerate in mass and the superposition 12tj|) is furthermore a solution of the mass equation JSJ. Indeed, as already 
noticed, S^tp^ is not an eigenstate of A 2 but a superposition of different A 2 eigenstates. On the other hand, the 
commutation relation [S^, P„] = implies S^tj)^ to have the same mass than <j)^ a \ This is possible only if the masses 
of different states are equal. 

Due to the Fock-space truncation, or to some other kind of approximation, the Poincare algebra is violated. The 
eigenstates (fM' are no longer degenerate and the solution l|26|) . built with eigenstates of different mass, cannot satisfy 
equation (JSJ). However, while this equation is an approximate one, the form (|26|l for the state vector remains valid. 
Each term in (|3 II) is an exact solution of the truncated mass equation 10 with eigenvalue M 2 . Their superposition 
satisfies no any mass equation but has the proper form of the non-truncated hamiltonian problem. The corresponding 
mass squared - at the same level of approximation - is given by: 

J 

M 2 =^c>4 2 (27) 

a=0 

The ensemble (0, M) obtained that way, constitutes the solution of the problem compatible with the degree of 
approximation considered. 

This formalism is translated to J = 1 states in the two-body sector as follows. The interaction kernel K(k' , fc, n, M) 
depends on scalar products of vectors k' ,k,n and also on scalar products with Pauli matrices: k-a, k ' -a, n-a. Therefore 
the total angular momentum operator constructed as 

J = —i[k x dj: ] — i[ft x d. n ] + s i + s 2 , (28) 

commutes with the kernel (s*i,2 are the fermion spin operators). In the cm. system this operator is proportional to 
the kinematical Pauli-Lubansky vector W given in (|23[l . The solutions of equation (JjJJ correspond to definite J and 
J z eigenvalues of the operators J 2 ',J Z . 

Since A 2 is applied to states with definite p, it has the form: 

A 2 = {n-jf. (29) 

A 2 commutes with the kernel K since J commutes with K and ft is a parameter. It commutes also with J since A 
is a scalar. Thus, like in the case of a full state vector lj2lfl) . the truncated solutions in the two-body sector are also 
labelled by a: 

A 2 if {a \k,n) = a 2 ^ a \k,n). (30) 

and the two-body wave function is a superposition of A 2 eigenstates ip^ with different a values: 

if(k,n) = c Tp (0) (k,ri) +ci ^{k.n). (31) 

The mass equations determining the eigenstates ip^ with different a are decoupled; in particular, the J = state 
is determined by one single equation. We would like to comment here that the decoupling into subsystems takes place 
in any formulation of LFD, both in the explicitly covariant and in the standard one. However, in the latter approach 
it looks as a matter of art, whereas in ECLFD this splitting has transparent reasons. For example, in |46| the four 
equations system for the wave function components with angular momentum projection m = was split, by a proper 
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transformation, in two subsystems with two equations each. In ECLFD this corresponds to the a = eigenstate of 
J = and J = 1 states, each of them having two components. 

Because the truncation of the Fock space, the states are not degenerate. Their splitting was effectively 
calculated in case of scalar particles in [44L \6nL l6fi| for J = 1,2 as a function of the coupling constant. It has been 
shown in |44| that this splitting indeed decreased when the interaction kernel incorporates larger number of particles in 
the intermediate states. However, the number of states taken into account in any practical calculation will be always 
very limited. The splitting, though decreased, will remain, specially for strongly bound systems like qq mesons. The 
problem of determining the state vector at a given level of approximation is thus not solved by this way. These are 
some of the reasons why, as explained before, our approach to deal with this problem follows a different philosophy. 
Despite the non degeneracy of ip( a \ we search the physical two-body wave function in the form l|3I|) . the same as for 
the full state vector l|2tj|) . The corresponding mass squared is given by: 

M 2 = CoA/ 2 + c\M 2 (32) 

where M 2 is the mass associated with ip( a \ The M 2 value thus obtained is always between Mq and M 2 , where would 
be the exact solution. 

To determine in practice coefficients c a , we use a method proposed in EH Iz2 > without explicitly solving J3UJ). 
It is based on the fact that, when the momentum k — ► 0, the interaction part in l|20(l is irrelevant and the angular 
condition reads simply L^ v (j) = 0. Thus, in this limit, ip does not depend on the light-front direction n anymore. Such 
a requirement unambiguously determines the coefficients of the superposition. The method was applied to a model 
with scalar particles |6a | and found to give very accurate results. The procedure will be detailed in section Ivl and 
illustrated by numerical calculations in section IvTTl 



IV. J = STATES 

The J = + two-fermion wave function can be written in the form [5?l l5^ | : 

&cr 2<Tl (ki, k 2l p,uJT) = ^/mu <T2 (k 2 ) <p U c u ai {kx), (33) 



where 



f 



a, 



,{k) = Ve^Tm g.k w ai (34) 



(s k + m) 



is the Dirac spinor normalized to u a u a ' = 2m5 aa ', w a the Pauli spinor normalized to wlw a > = 8 aa > and Ek — 

+ m 2 . U c — 7 2 7° is the charge conjugation matrix. In its turn, <p is written as a superposition of two independent 
spin structures Si 

s > - ik:- (35) 

£ k ( 2m A m 2 \ 

<J2 = — ?= u 9- 75 

2V2mk sin 9 \u>-p e% J 

whose coefficients /i, scalar functions depending on variables (k,cos9 = n-k/k), are the wave function components in 
the spin-space: 

- fiSi + f 2 S 2 , (36) 

The existence of one additional component with respect to the non-relativistic theory is due to the cD = uj^j^ term. 
The number of independent amplitudes determining the wave function is however the same, whatever will be the LFD 
version used. We have shown in a preceding work [68( that the J w = + state we are considering is strictly equivalent 
in the standard approach to the (1+, 2—) one ^(| which is described also by two components <E> 1+ , <I> 2- . 

In the reference system where k\ + k 2 = the wave function ilM-il) takes the form: 

$(72(7! = Vmwl 2 ip(k,n)wl 1 (37) 
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with 

The definition of the components themselves is to some extent arbitrary, as are arbitrary the choices of structures 
(I35J1 . Our choice (|35[1 is justified by the clear separation of n-independent and dependent terms it induces in the wave 
function jggjl. 

The normalization condition reads: 



1 f , ^ ,2 d 3 k m /" r -, f , , , f .,d 3 k 



CT2CT1 



(27r) d 7 s k (27r) J J e k 

where we denote (f> = 7o<^7o- The spin structures Si introduced in (|35|l are orthonormalized relative to the trace: 

Tr{S i (k 2 +m)S j (k 1 -m)}=6 ij , (40) 

where Si = joSj^fo, that is: 

S = £k / 2mc2> m 2 \ 

2 2V2mfcsin6< 75 \ w-p e\ ) 

Substituting in (J5J the wave function (|33|l . multiplying it on left by u(k 2 ), on right by u(k\) and using relation 
J2 a U<J (k)u a (k) = k + m, we find: 

[4(£ 2 +TO 2 ) - M 2 ] (fc 2 + m)0(fci - m) (42) 
to 2 /" 1 ,« , f , - ,d 3 /c' 



27T 3 7 Am 2 {Q 2 + ^ + m )°^ k 2 + rnW(k[- m )0 1 (k 1 - m )- £ki 

with O = U c O t U c - Replacing here <j> by its decomposition (|36|1 . multiplying equation l|42|) by ft and using the 
orthogonality relations (|40|l . we end up with a two-dimensional integral equations system for components 

[4(P + m 2 )-M 2 ] /«(M) =-^Lj2f K^kAk'^mk'^') — , (43) 

Its solution will directly provide the mass of the J n = + state. 

Kernels appearing in (|43|) result from integrating kernels k^ over the azimuthal angle ip: 



Kij dip' 
m 2 e k e k > J Q 2 + 2ir ' 



K ■ - t / —221 2tL (MS 



with Q 2 defined in (|12|) . For S, PS and PV couplings Oi = U c O[U c = 0\ and fty are given by 

Kij = \e k e' k Tr\s i {k 2 +m)0 2 {k l 2+m)S'Ak\-m)0 1 {k l -m)\ (45) 
4 L 

We denote by S'j the quantities (|35|l as a function of primed arguments. For vector exchange: 

Kij = --e k e' k L aP Tr \s,{k 2 + m)0%(k' 2 + m)S' j (k' 1 - m)0((k ~ m)j , (46) 



Tensor L a p is defined in l|14l) and we have taken into account that for V coupling 0\ — U c O\U c = —0\. The analytic 
expressions of Kij for S, PS, PV and V exchanges are given in appendix lAl 

One would remark that we have kept, for convenience, a three-dimensional volume element in equation (|43|l despite 
the fact that kernels as well as amplitudes fj are independent of variable ip' . 
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V. J = 1 STATES 

In a similar way than in (|33[) . the J = 1 + two-fermion wave function can be written in the form |56ll58|]: 

^^(fci^.p.wr) = e M (p,A) S CT2 (fc 2 ) V U<M^{k x ) , (47) 
where e„(p, A) is the polarization vector. <^ develops over the six spin structures 

c _ (^1 - k 2 r 1 _ (fci - * 2 )"6 

— n 5 ) ^2fi — — 7 i — , ^4^ — ~ ) 

Zm z m w-p Zmuj-p 

S 6lt = ^75e Ml/p7 fci^ 2pW7 , 5 6AI = ^£ (48) 

with components (pi, invariant functions depending on the same scalar variables than for J = 0, 

4>^ = ipxSift + (P2S211 + <^3<5 , 3m + ViSifi + psS^ + (feSeti ■ (49) 
In the reference system k\ + k% = this wave function takes the form: 

$„ 2ai (k,n) = y/mwl 2 ^(k,n)a y wl x , (50) 

with 

$(k,n) = h ^ + f 2 \{^^-a) + h 1 -{?,n{n-S)-a) 

+ f^(3k(n-a) + M{k-a) - 2{tn)a) + h\[\\[k x n] + f 6 ^(H{k-a) - k(n-a)) . (51) 

Contrary to the J = case, components appearing in l|51[) are not the same than ip t from (|49|l . Their relation 
is given in Appendix iBl Components fa 4,5,6) driving n-dependent spin structures, are of relativistic origin and are 
absent in a non relativistic approach. 

As explained in section ITTT1 the system of equations determining the six components fi is split in two subsystems, 
corresponding to the eigenvalues a = 0, 1 of A 2 (|29|) . Like for the J = wave function, the J = 1, a = eigenstate 
is determined by two components whereas the remaining four correspond to J = l,a = 1. We would like to notice 
that the total number of components as well as the dimension of decoupled subsystems (2+4) coincides with what is 
found in the standard approach |4fij . 

The components determining the eigenstates of A 2 will be respectively denote by fl , }=7 < 2 an d 3i=7 2 3 4- They are 
indeed different from the fi appearing in the wave function (|51|l though g's fully determine /'s by linear combinations. 
In view of constructing the superposition (|31[1 it is convenient to represent the eigenfunctions m the form l|51|l . 
Only some of the six f^ involved components will be independent - two for the a = state and four for a = 1 - but 
this way of doing will facilitate further analysis. 

In the two coming subsections we will explicitly construct the eigenfunctions ip^ of the kinematical operator A 2 , 
obtain the corresponding mass equation (JSJ) in terms of g\ a ^ and relate them with components f^ defined in l|51|l . 



A. a = 

One can check from equation (|3U[1 that tp^ is parallel to ft, i.e. it satisfies ip^ = nin-ip^), and has the following 
general decomposition: 

It can be written in the form (|51|l by defining the components 



1 - cos d 9\ - ~y= sm <?92 
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A 0) = o 

f (o) _ 



AO) _ 
Ji — 



1 



V2 
%/3sin0 



(0) 



(0) 



fh 
1 



cot ( 



,(0) 



p(0) 
5 

(0) 



(53) 



that is four non-zero components, with only two of them being independent. It can also be represented in a four- 
dimensional form similar to (|49|l 



by introducing the spin structures 



CO). 



AO) c(o) ,(o) q (0) 



(54) 



5' 



(n) 



o(0) 



with S ifl , defined in (gSJ. 

The normalization condition is 



V3Af 
2V2fc 3M ' 
V3Af /to 2 cos 



i\/2sin6) V 2e fc fc 



5V + I ! 



3(2tt) 



3(2tt)3 



— V /V I 

!tt) 3 ^ y 1 ct2<ti1 



d 3 fc 



£fc (2tt) 3 



W u Tr{cj)W(k 2 + m)^ 0) (fc a - to)} 



Tr{^ (0)t (k,fi)^(k,fi)}— = 



(0), 



d 3 fc 



£fe (2tt)3 



(. 9 i 0) ) 2 + ( 5 f) 2 



d 3 k 



= 1. 



with 



(55) 
(56) 



(57) 



(58) 



The spin structures S^J are orthonormalized relative to the trace operation in (|57|) (cf. eq. (|40I0 : 



n^Tr{^ } (fc 2 + m)S^>(k 1 - to)} = 5 
Note that <Sj^ = 7 S' l ( ° )t 7o = S$ . Similarly to equation g2J we get: 



(0), 



(59) 



4(fc 2 + to 2 ) - M 2 } (k 2 + to)4 0) (fci - to) 



2tt 3 y 4to 2 (Q 2 + ^ 2 ) 



(fca + m)0 2 (fc 2 + ™R (0) (*i - m)6i (fci - to) 



ri 3 fc' 



(60) 



In order to obtain the system of equations for components , we multiply equation 1)60(1 by IP 1 ' and . Taking 
the trace and using the orthogonality condition 1)59(1 we obtain the system of equations: 



4(P + to 2 ) - M 2 ] gf (k,n) = -|L J %%' \fi)gf\k' ,n)^£ 



d 3 k> 



(61) 



i=i 



which provides the mass of the J = 1, o = state. They have the same form than 143|) . with kernels JQ- given in 
terms of k-^ integrated over the azimuthal angle ip': 



K ia) = — 

*3 ™2 



1 



2tT JO) , , 

K ij dip' 
> 2 + n 2 2tt ' 



(62) 
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For S, PS and PV couplings they read 



4 0) = \e k e' k W v Tr \s^(k 2 + m)0 2 (k' 2 + m) S'}? (fc'i - m) Oi(fei - m) 



(63) 



where denotes (|55l) with primed arguments. For vector exchange 



4j> = --e k e' k U^L af3 Tr [S^(k 2 + m)0%(k' 2 + m)SjJ(k[ - m)Of(*i - m) 



(64) 



Tensors L a/ 3 and IP^ are defined in (|14fl and l|58|l . The analytic expressions of nf) for S, PS, PV and V (/ t = 0) 
exchanges are given in appendix lAl 



B. a = 1 

It follows also from i|30[) that the A 2 eigenfunction corresponding to a = 1, is orthogonal to n, i.e. satisfies 
n-t/jW = 0. To fulfill this condition, it is convenient to introduce two vectors <7j_) orthogonal to n: 

? k — cos 9ft 

fc -L = : — a 

sintf^ 

a± = a — (n-(f)n. 



with fc = k/k and | kj_ | = 1. Function i/A 1 ) obtains then the decomposition, analogous of H52|) . 
^ (1) (fc,n)=. 9 ( 



(65) 



in terms of the four scalar amplitudes g^- It can also be represented in the form (|51JI by defining components /, 



, AD 



/i 1 
/I 1 



/i 1 



/ 5 {1 



2 (i) 
2 



-o (1) 



\/3sin 2 (9 

1 (1) (1+ COS 2 6>) (!) x/2 



2\/3( 



\/3sin 2 



-52 



V3 



cot( 



o ■ 2 a "52 

3sm # 



(i) 



\/6sm6> 



o (1) 



sin0 



.94 



(1) 



and in the four-dimensional form 0^ similar to (|50|l : 



_(l) c(l) , „(1) c(l) , _(l) c(l) , _(l) c(l) 
9i •Jifj, + g 2 d 2fi + 53 a 3n + 9a ^Ah ■ 



The four spin structures S^J are orthonormalized according to (|59|) and read: 



S il = h v S i^ i = 1, . . . , 4; j = 1, . . . , 6, 



(66) 



(67) 



(68) 



with Sj M defined in (|48|l and /ly coefficients given in appendix [5] The normalization condition in terms of cp^/ and 



ip^ exactly coincides with (|57|l . In terms of components g^ it becomes: 



(5l 1) ) 2 + (^ 1) ) 2 + (^ ) ) 2 + (5i 1) ) 2 



(2tt) 3 



d 3 k 

£k 



(69) 
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The system of equations for the scalar functions g\ is obtained similarly to l|61|> and reads 
4(P + m 2 )-M 2 

It is the mass equation of the J — \,a — 1 states. Kernels K-j are calculated in a way similar than 1)62(1 . The 

corresponding are obtained with the replacement — > in l(63|l and (|64|) . Their analytic expressions for S 
and PS exchanges are given in appendix 1X1 



g f\k,fi) = -^ ^K^\lk',n) 9 ^(k',H)^-. (70) 



C. Physical solution 

The solutions tp^ constructed in the preceding sections, although being exact eigenstates of the truncated Hamil- 
tonian, are only auxiliary. As explained in section ITTT1 the solution satisfying the angular condition (|20H is given by 
the superposition 1(3 1|) of states with different a. The coefficients c a of the superposition can be obtained by solving 
the angular condition in the truncated Fock space. We will show in what follows that they can alternatively be 
determined by imposing the independence of the wave function on the light-front vector ft at k = 0. 

In order to do that, it is convenient to write down ip^ in the form 1(51(1 with the components /5 given by equations 
((5311 and l(66|) . Written in term of /'s the superposition l|31|) reads 

/i = co/i° ) +c 1 /P. (71) 

The condition that ip(k = 0,n) does not depend on ft becomes: 

dgfi{k = Q,e) = i = 1,2 (72) 
/,(fc = O,0) = j = 3, 4, 5, 6. (73) 

Let us show that there exists two coefficients c a , normalized to c§ + c\ — 1, satisfying the above six equations. They 
are determined by the only values at k = of the first components . 



To this aim, we consider the behavior of f- (k,z) in the k — > limit. The components in front of structures 

.(a) 

2,4,5,6- 

f2A, 5 , 6 ( k = > d ) = a = > 1 - ( 74 ) 



involving the unit vector k are f^l 5 6 - By construction, they must vanish at k = 0, i.e. satisfy 



Concerning a — states, this condition is trivially satisfied by /| 5 since from Ij53(l . they are identically zero whereas 



12,5 

f$ will satisfy O if: 

^ 0) (A; = 0,e) = +6 cos6», (75) 
^ 0) (fc = 0,6») = -b sin 9, 

bo being a priori an arbitrary function of 9 which later on will be shown to be constant. The only components which 
are non zero at k = are f[ ^ . Inserting the values l|75|) in l|53|l we find: 

A (0) ((M) = ± COS 9g?\0,9) - ±sm9gi 0) (0,9) = -j=b 



Concerning a = 1 solutions, determined by four independent components g^\ we see from (|t)6^ that condition l|74|) 
implies 5^ 4 (0, 9) = 0. The only non-vanishing component at k = is thus 5^ and we will denote by &i its value: 

g{ 1 \0,6)=b 1 . (76) 
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By inserting this values in (|66() we get: 

/i (1) (o,0) - vf 5 * 1)(M) = vf 61 ' 

/ 3 (1) M = -^(O.g) + ^+g^ ^(0,g) - ^cot^(O^) = --L fcl . 

Components /g ^ and /g 1 ^ are the only n-dependent structures which gives non-zero contributions at k = in the 
corresponding wave functions ip^ and if)' 1 '. These contributions must cancel in the physical wave function ip, what 
gives the relation 

/ 3 (0, 9) = co /f(0, 9) + Cl / 3 W (0, 9) = co^bo - c^h = 

This relation, together with the normalization condition c\ + c\ = 1 , allows us to determine the coefficients c a of the 
superposition l|31(l . They read: 

We see from the above expressions that conditions (|72|l and 1)73(1 will be satisfied if and only if coefficients b a are 
actually independent of 9. 

It is worth noticing that if the wave function ip does not depend on n, these coefficients becomes especially simple: 

co = yf, C1 = . (rs) 

Indeed, from an n-independent wave function -0 we can construct normalized n-dependent states with definite a as 
follows: 

tf^ \k,n) = VZn(n-$(k)), 

The initial function ij) (k) is reproduced by taking their superposition with coefficients (|78H . In the case of scalar 
constituents, we found |63 that coefficients c a are very close (with the accuracy w 1%) to the values J7SJ), despite the 
fact that the wave function strongly depended on n and the split between Mq and Mi masses was large. 

Let us finally summarize the procedure followed to construct the physical wave function. The solution of the mass 
equations l|61|l and J7DJ), provides the mass squared M% and the components gfl and gi_4 of the A 2 eigenstates. The 

non-zero values of the first-components g[ a ^ at fc = determine - by means of (|75|l and (|76l) - the coefficients b a . 
These are inserted in equation l|77l) to provide c a , coefficients of the linear combination determining the physical mass 
M 2 from (|32|l and the components (|71() of the wave function l|51() . Components of this superposition are related 
to by lf55T) and correspondingly. 



VI. NONRELATIVISTIC LIMIT 



In the forthcoming sections the LFD results will be compared to the corresponding non relativistic limits. We mean 
by that the zero order terms in the 1/m expansion of the LFD equations and kernels. This section is devoted to 
precise how this limit is obtained in the different OBE kernels, having in mind in each case (i) what are the LFD 
wave function components that should be retained and (ii) what kind of equations will they satisfy. 

In order to have some insight in the weak coupling limit, but also as a test for numerical calculations, it is often 
useful to consider the LFD solutions as a perturbation of the non relativistic wave functions. This approximation 
was used in [5^, |5|j to calculate the NN S-wave function and deuteron electromagnetic form factors [6lJ . We will also 
present in what follows how these first order relativistic corrections can be obtained in the different mass equations 
Pfy. (E2) we consider. 
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A. J = states 

For the scalar exchange the leading contribution in the kernel matrix is, according to (|A4|) : 

Kn = --H^T 2 = V s (k-k>) (79) 

(k - fc') + i± 2 

Corrections to this kernel are of the 1/m 2 order both in diagonal and non-diagonal terms. It follows that the J = 
wave function lj38[) . contains in the non-relativistic limit the f\ component only, which is furthermore independent of 
0. Introducing non-relativistic kinematics, i.e. 4(fc 2 + m 2 ) — M 2 w A(k 2 + mB) where B — 2m — M << m is the 
binding energy, the equation for fi = Jnr component becomes: 

d 3 k' 
(2tt) 3 



(k 2 + mB)f NR {k) = -m / V s (k- fc')/iVfl(fc')77TZW ( 80 ) 



with kernel (|79|l . This is the Schrodinger equation with the Yukawa potential Vs(r) = — aexp(— \xr)jr. 

For vector exchange we obtain the same equation H8()(l with a kernel differing from 1|79[) by a global sign. This 
corresponds to the repulsion between two fermions (e _ e _ , for instance). 

We see that for the scalar and vector couplings, the non relativistic limit of LFD equations coincides with the 
one-component Schrodinger equation. 

For pseudoscalar and pseudovector exchanges the leading diagonal kernels are of the 1/m 2 order, whereas the non- 
diagonal ones are of 1/m 3 . Thus, for these couplings the non relativistic limit does not exist. In the leading order 
and since the K22 kernel is repulsive, only the /1 component remains. The corrections due to are expected to be 
bigger than for scalar and vector cases. Component /1 satisfies at this order the Schrodinger equation l|8U|) with a 
kernel proportional to 1/m 2 : 



v PS (k-k') = — = ™ 



{k-k') 2 +ii 2 
In coordinate space it corresponds to 



M 2 



(k - k') 2 + fi- 



/j, exp(— fir) 
Ait r 



(81) 



(82) 



For these couplings the leading term is of the same order as relativistic correction in the scalar and vector cases. We 
will see that a similar situation takes place for the J = 1 state. This fact makes an important difference between the 
couplings. Pseudoscalar and pseudovector exchanges appear always as being relativistic corrections. 

We would like to remark from the above results that in the non relativistic limit the n-dependent terms in the LFD 
wave function (|38|l and kernels disappear. 

For models involving the sum of all exchanges (like for the OBE NN interaction) the non-relativistic limit is 
determined only by the S and V exchanges. First order corrections can be obtained by inserting the non-relativistic 
component fi — /nr into the r.h.-side of equations l|43() . 

m 2 P ii 3 k' 
[4(k 2 + m 2 ) - M 2 ] fi{k, 6) = -— / K a (k, 0; k' , 6')f NR {k') — . (83) 

They generate a perturbative solution for the two components which incorporates the first order relativistic effects. 
This approach was followed in [5^ to obtain the 1 Sq NN scattering wave function. 



B. J = 1 states 

For J — 1 states, components g\ obtained by solving the mass equations differ from those appearing in the wave 
function (fi). Our first step is to determine the form of g\°^ in case of a non relativistic wave function. The non 
relativistic wave function components do not depend on n and, according to (|78[) . are given by: 
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Substituting (|84|) into (|53|l and (|66(l we obtain a relation between fi and g" 1 components. These equations are solved 
relative to g^ and the result, expressed through fi, reads: 



A 0) 



9 f ] 



9 P 



9 P 



(1) 



/i cos 9 + /2V2 cos 9 + /aV^cos 6 + f< 



7 + cos 29 
2V2 



+ M 7; sin 



-fx sine + f 2 ^= sin9- / 3 V2sin0 - / 4 ^= sin 29 + f 6 ^, 



1 29 



l + 3cos2fl 1 r- 

/l_/a— Ivi /3 71 _/4 ^ 



cosf 



= /a 



2\/2 



sin 



3 3 \/3 

/ 2 - sin26» + fi- sind- / 6 — sin( 



/5\/ 2 sin( 



(85) 



As previously discussed, in the non relativistic limit there are no n-dependent terms in the LFD wave function l|51(l 
and only fx and f 2 components among the six fi survive. We have shown in that one actually has f\ ~ us, fi ~ 
—lie /3-6 ~ where us and ud are respectively the usual S- and D-wave non relativistic components. Inserting their 
expressions in (|85|l we obtain the form of the non-relativistic functions g: 



9? 



A 1} 



9 P 

9? 



(us — V2ud) cos 9, 
1 



US + UD I SID fi. 



1 



u s 



4V2 



(1 + 3cos26>) u D , 



2V2 



sm v ud, 



- — sin 29 ud, 



= 0. 



(86) 



We see here that the ^-dependence of the auxiliary components g!- remains even in the non relativistic limit. It will 
disappear only in the linear combination giving the physical components f\,i- 

Let us first consider the scalar exchange. The mass equation for a = eigenstate (|61|) and the scalar kernels IjA If) , 
becomes in the leading order (l/m)°: 



C{k)gf\k,9) 



-Aan cos 9 



g^ (k', 9') cos 9' - g ( z> (k' , 9') sin 9' {. . .} 



C{k)gf\k,e) = +4a7rsin6» / \gf\k', 9') cos9' - ^ 0) (fc', 9') sinfl'l {. . .} 



d 3 k' 

£k> 
d 3 k' 
£k> 



(87) 



For shortness we denote by C(k) the kinematical part and by {. . .} the kernel contributions which are common to all 
couplings and states. 



{...} 



1 



1 



m 2 e k e k > Q 2 + fi 2 



These factors contain 1/m and 1/m 2 terms but we do not write them explicitly and analyze only the kernel contri- 
butions resulting from Kij . 

Since the integrals in the right hand sides of (|87J) are the same, its solution has the form: 

g { °\k,9) = +gW(k) cos 9, 
gi°\k,e) = -g^(k) sine. 
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with g(°^ an unknown function to determine. Substituting l|88|) into 187(1 we find the equation for g^: 

/H 3 k' 
gW(k'){. ..} — (89) 
Efc' 



For a = 1 state, we found in a similar way that only g^ survives and satisfies to the same (l/m)° order, the equation 

C(k)g?\k,9) = -4cm f ff) {. . .} — (90) 



It coincides with the equation l(89() for a — and, hence, provides the same mass. We see in this way that, in the 
leading order, a = and a = 1 states arc degenerate. The coefficients c a of the superposition 1(71(1 are calculated in 

terms of b a given by 1(75(1 and ((76(1 . Since in the leading order, g^ and equals <?(°), one has bo = b\ = g^(0) and, 
from 1(77)1. the values f78)). 

In next to leading order - 1/m 2 - we get for a = 0: 



C{k)g{ 0} = air-^ [ \2g^kk' — 3(k 2 + k' 2 ) cos9(g'^ cos9' — g'S ^ sin#')l {. . .} 
m z J l ' i 



d 3 k' 



C(k)g { 2 0) = mrsin^ J 3(fc 2 + k 12 )^" cosO' - g'^> sin 5') {. . .} (91) 
and for a = 1: 

C{k)g^ ] = J [3(k 2 + k' 2 )gf ] + kk' cos9{-2g'^ cosff + V2g' 3 W sin 00] {. . .} ^ 



£k> 

d 3 k' 



C(k)g^ ] = -mrsin(9^- I kk'(V2g\ K1 > cos9' - g' {1 > sin9') {. . .} (92) 
'"- ./ ' " £k' 

These systems of equations - ((9111 and ((92(1 - are already different and the masses of the two eigenstates are split. 

For vector exchange, the situation is quite similar. The equations in the leading order (1/m) differ from 1(87(1 an 
((90(1 only by a global sign in their right hand sides. Thus for these two couplings, as it was the case for J = 0, the 
leading order is (1/m) . 

For pseudoscalar exchange, the leading contribution in the kernel has order 1/m 2 . Indeed, from the analytic 
expressions given in I A 21 we found for the a = state: 



3in e-L f kk'iV^g'^ cos 9' - g 3 {1) sin 9') {...} — 
m 2 J ' ■' s k > 



(0) 

9i 



air cos 9-^ J k 12 (gf ] cos 9' + sinfl') {. . .} 
gf ] = -mrsinfli J k' 2 (g^ cos 9' + gf ] sin 9') {.. .}^- 
Like for the scalar coupling, the solution of ((93(1 has the form (188(1 with g^ satisfying the equation: 



d 3 k' 
£k> 
d 3 k' 



(93) 



(o) = an JL f fc ' 2 ( cos 2 ff _ sin 2 e > )g (o) { _ } ^ (94) 



(95) 



m J £k 
For a = 1 the leading order equation reads: 

d 3 k 



^=-^/*'W^f>{...}^ 



which is now different from 194(1 . The masses M and Mi calculated with pseudoscalar exchange are therefore always 
different. Their difference remains even in systems having small binding energies or when the large momentum 
contributions are removed using small cutoff parameter A in form factors ((18(1 . 

Pseudovector exchange kernel differs from the pseudoscalar one by the replacement 75 — > 75 — "J 57 " or by 75 — > 

75 + ( see ec l - (4-18) m There is so an extra term proportional to "Im oc fc which does not contain 

(l/m)° terms. The situation is therefore the same as for the pseudoscalar case. 



18 



To summarize, we have shown analytically that in the non-relativistic limit for scalar and vector exchanges, the 
energies B(a — 0) and B(a — 1) coincide with each other and the coefficients cq, c\ tend to y| respectively. On 
the contrary, for pseudoscalar and pseudovector couplings this is not the case. In this sense, for the pseudoscalar and 
pseudovector exchanges, the non-relativistic limit does not exist. If the kernel is the sum of all the exchanges, like 
NN kernel, the situation is the same as for the scalar and vector exchanges, since in non-relativistic limit the (l/m)° 
order dominates, resulting from these exchanges. The existence of deuteron, for example, as a nonrelativistic system 
(with a reasonable accuracy) is due to contribution of the scalar and vector exchanges in NN interaction. 

Perturbative solutions are obtained by substituting the zero-order functions i|86|) into the right hand sides of LFD 
equations l|61ll and (|70() . If D-wave is neglected, the six perturbative components are given in terms of the only non 
relativistic wave function us simply by: 



4(fe 2 +m 2 )-M 2 l .9 4 (0) (fc,n) 
4(P +to 2 ) - M 2 ] ,9 4 (1) (fc, n) 



m 

'27 

—J / K ii u s{k ) 



d 3 k' 

£k> 



2tt 3 



Efc' 



(96) 
(97) 



We would like to mention here that one appreciable advantage of the LFD formalism with respect to other relativistic 
approaches is the clear link it has with the non relativistic dynamics. On one hand because LFD wave functions have 
the same physical meaning of probability amplitudes. On the other hand, because their components fi split in two 
families: those which in the non relativistic limits become negligible and those which tend to the usual non relativistic 
wave functions. 



Next sections are devoted to show the numerical solutions obtained with diffcrents couplings. Their very different 
behaviour motivates to be treated separately. 



VII. RESULTS FOR SCALAR COUPLING 

Our first results concerning the Yukawa model have been reported in |68l 169^ . The main interest in these papers 
concerned the stability of the J = 0,1 solutions with respect to the cut-off, i.e. the possibility of getting stable 
results without any vertex form factor. We showed in particular that J = + states were stable for coupling constant 
smaller than some critical value a < a c ~ 3.72 and unstable above. On the contrary the J = 1 + states were found 
to be unstable for any value of the coupling constant and both projections a = 0, 1. This instability manifests in the 
logarithmic decrease of M 2 (k max ) for a given value of a - or equivalently of a(k max ) for a given value of M - and 
imposes the use of form factors. 

We first consider the J = + state. Its wave function is determined by two components /». Although the use of 
vertex form factors (FF) is not required [68|. we would like to notice that the convergence as a function of k max is 
very slow. Unless otherwise specified the results that follow correspond to /i = 0.15. 

For a weakly bound system (B=0.001), the coupling constant found solving LFD equations is ay tifc=0.331 whereas 
the non relativistic (NR) value is ajvi?=0.323. By the latter we understand, the results obtained by inserting into 
the Schrodinger equation (|80|l the static potential (|79|l resulting from the leading order approximation as has been 
discussed in section IvTl Like in the Wick-Cutkosky (WC) model - scalar particles interacting by scalar exchange - 
relativistic effects are repulsive j64[. They account for only a 3% difference in the coupling constants whereas in WC 
they are sizeably bigger (avKC=0.364). 

Corresponding wave functions are displayed in Figs. QJl and0J}. One can see that component /i dominates over 
/2 in all the interesting momentum range and that fi has a zero at k ss 0.25. One also notices in Fig. that f\ is 
very close to the NR wave function in the small momentum but it sensibly departs with increasing k; for k ~ 1.5 the 
differences represents more than one order of magnitude in the probability densities. The coupling between the two 
relativistic amplitudes has a very small (0.1%) attractive effect in the binding energy. 

In the strong binding limit (B=0.5), the situation is quite similar with enhanced relativistic effects in binding 
energies and wave functions. One has ay u/ t=2.44 for unr=1.71 and the differences in the wave functions - displayed 
in Figs. and|Sj3 - are already visible at k — momentum (Fig. One can see however in Fig. that - even for 
deeply bound systems - f± component still dominates over fa. 

It has some interest to compare the LFD results for Yukawa (two-fermion) and WC (two-scalar) models with the NR 
results. We have displayed in Fig. [S]the corresponding coupling constants for different values of the binding energy. 
One can see that the Yukawa results (ayuk) are systematically closer to the non relativistic values than awe are, as 
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FIG. 5: LFD wave function components fi for scalar coupling (B=0.5, n = 0.15) in linear (a) and logarithmic (b) scale compared 
with the non relativistic solutions 



if the fermionic character of the constituents generates closer binding energies to the NR ones but larger differences 
in the high momentum components of the wave function, due to the different asymptotic of interaction kernels. 

Though not necessary to get stable solutions, form factors they have been widely used in most of the preceding 
OBEP calculations performed in momentum space |60j| . It is thus interesting to estimate their influence in the 
predictions. To this aim we have considered the vertex form factors used in the Bonn model (|18f) with, for the 
scalar coupling, n=l and A=2.0. Their effects are found to be repulsive. For B=0.001 they remain relatively small 
(aYuk— 0.376 instead of ay u fc=0.331) but for B=0.5 the differences reach already a factor two (a^ u fe=5.32 instead of 
otYuk = 2.44). It is worth emphasizing that whatever will be the degree of refinement in the dynamics, the results of 
a relativistic calculation will be strongly influenced by this phenomenological and not well controlled trick. 

The system of equations for determining the J = 1 + , a = (%=72 ) anc ^ a ~ 1 G?i=7 2 3 4) solutions are both unstable 
and require cutoff regularization |4rll68|. This can be seen in Fig. where the ct(k ma x) variation for a = and a = 1 
cases displays a logarithmic dependence. One can also see in this figure the non degeneracy of both states due to 
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FIG. 6: Comparison of B(a) between the Yukawa (dashed line) and Wick-Cutkosky (dot dashed lines) models in LFD and non 
relativistic (solid line) results in J = + state 



the Fock space truncation discussed in section IIIII We remark however that if the binding energies - or equivalently 
coupling constants - of states with different projections a are not equal, they are almost-degenerated in a wide range 
of k max values. For instance, at k max = 10 one has a a= o=1.17 and a a= i=1.18 while at k max = 90 one has a a= Q=1.14 
and a a= i = 1.16. These weak splitting - of less than 1% - for a noticeably bound system (B = 0.05), are rather 
surprising in view of the results obtained in the purely scalar WC case |65l 66] , in which the difference in coupling 
constants for the same binding energy is rather 20%, what corresponds to AB w B. 

The <7. solutions for a=0 and a=l states are respectively represented in Figs. and [§] for several values of 
0. They were obtained with a coupling constant a s = 1.18 and a sharp cutoff at k max = 10. We remark that 
with the conventions used g% (k,0) — and on has g^(0,0) = — g% (0, 90°), as expected from (|75f) . In addition: 

<7i (0,0) = — 52°'(0, 90°) ~ gi^iO, 0), as expected from <|76[) and from the fact that coefficient Cj, defined in (|77|l . are 
very close to the values <|78|) . Corresponding binding energies are B a= o — 0.0506 and B a= \ = 0.0498, values which 
are 1% close to each other. The splitting of the binding energies is an increasing function of the coupling constant. 
Figure ITUl shows the calculated B a (a) dependence for both J=l eigenstates. For a s = 0.55 the values are respectively 
B a=0 = 9.7 10" 3 and B a=1 = 9.6 10" 3 whereas for a s = 2.87 B a=a = 0.523 and B a=1 = 0.467. The non degeneracy 
remains reasonably small even for strongly bound systems. 

The six components fi of the J = 1 + physical wave function are determined by a linear combination l|71l) of 
functions which in their turn are expressed in terms of g\°^ by (|53|l and (|66|l . We remind that coefficients c a of 

this linear combination are computed from components t^ 0-0,1 ^ only. For the solutions presented in Figs. [H] and O 
they are found to be Co = 0.582 and c\ — 0.813 and the corresponding energy is B=0.0501. Note that these values are 

very close to those obtained in case of rt-independent interactions (|78|l : c o — ^ — 0.577 and c± = = 0.816. They 
become even closer to these values for smaller binding energies and they smoothly depart for strongly bound systems. 
For a state with B w 0.5 and the same sharp cutoff k max = 10 one has for instance cq — 0.610 and Co = 0.793. 
Components fi thus obtained are displayed in Fig. ^2f° r 9 — 30 degrees in linear (a) and logarithmic (b) scales. One 
can see that component /i dominates over all remaining five in all the momentum range. Among the components 
of relativistic origin there is not a clear dominance. Notice the very small value of fi component, corresponding to 
the tensor D-wave, that would be absent in a non relativistic approach. These components have a definite parity in 
variable cos#, /i,2,3,5 being even and f^g odd, as shown in Fig. Illb for a fixed value k = 1. 
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FIG. 7: Logarithmic dependence of the coupling constant as a function of cutoff for the J = 1 + a=0 and a=l states. Calculations 
correspond to B=0.05 and fi — 0.25 
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FIG. 8: <?f = i 2 solutions for scalar coupling with a = 1.18, /i = 0.25 and sharp cutoff at k max = 10. Binding energy is 
B = 0.0506. 



VIII. RESULTS FOR PSEUDOSCALAR COUPLING 



F or p seudoscalar coupling, the stability analysis was performed using the same methods than for the scalar one 
[7l|,[72| and presents some peculiarities. 

Equations for J=0 + states are found to be stable without any regularization. The asymptotic behavior of the 
pseudoscalar kernel is the same than the scalar one it has a repulsive character which do not generates instability. 
The results lead to a quasidegeneracy of the coupling constants for binding energies which vary over all the physical 
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range [0, 2m]. One gets for instance, a = 55.4 for B = 0.001 whereas a = 58.5 for a binding energy 500 times bigger 
B = 0.5, showing an extreme sensibility of this model to small variations of the coupling constant. The origin of 
this behavior was found to lie in the second channel equation (K22) and has been understood analytically 01 with 
a simple model. The use of form factors - though not required for the convergence of solutions - is necessary if one 
wishes to eliminate this unusual a{B) dependence. Calculations have thus been performed using form factors i|18|) 
with n=l and A=1.3 as in the Bonn model. 

In the weak binding limit (B=0.001) one has (Xlfd = 190 and cxnr = 166, a repulsive effect much stronger (15%) 
than in the scalar coupling. Corresponding wave functions are shown in Fig. 1121 One can see that the component of 
relativistic origin fa w f± at k ~ 0.3 and dominates above k=l. A similar result was found in the np 1 So scattering 
wave function calculated perturbatively with all the OBEP kernel in [5!j. Contrary to the Yukawa model, the role 
of relativistic components is crucial already for such a loosely bound system. The coupling between components is 
also very important: by switching off the non diagonal kernels K12 — K21 — the coupling constant moves from 
O-lfd — 190 to otLFD = 251. It has thus an attractive effect which tends to minimize the difference between LFD 
and NR results. The comparison between /i and the non relativistic solution Jnr shows a very good agreement in 
the small k. When k increases, large differences appear and Jnr has even an additional zero at k = 1.1. 

It is worth noticing the dramatic influence of the form factor in all these calculations. One has for instance 
(xlfd~W3 for A = 5 and aLFD — ^T25 for A=0.3. We remind that the value used in the Bonn model for this coupling 
IS A Bonn = 1.3. 

Quite surprisingly, in the strong binding limit (B=0.5) we have found a£i?£)=1462 and 0^^=3065. Relativistic 
effects become now strongly attractive (cxlfd < &nr)- An essential part of this attraction is due to the coupling 
of the two fx — / 2 components in the LFD wave function. By performing one channel calculations, one has indeed 
C(lfd=3001, what represents a strong reduction in the effect though it remains slightly attractive. We have checked if 
this attractive effect happens for different values of the exchange mass \x. For the same binding energy {B = 0.5) and 
fi = 0.5 we have found o.lfd = 1728 and aNR=^00, repulsive once again. It is worth noticing that for this coupling 
cxnr is a decreasing function of [i whereas olfd increases, at least in this energy region. This tells us the difficulty 
of talking about the "sign of relativistic effects" in general. They turn to depend not only on the kind of coupling but 
also on the binding energy of the system and - furthermore - on the mass of the exchanged particle. 

It is interesting to study the zero binding limit of the LFD results and compare them with the non relativistic ones. 
The NR potential (|82|l has been modified by including the Bonn form factor (|18H . The results are given in Fig. ITSl 
for an exchange mass [i = 0.5 and with two different cutoff parameters A in the form factors. They show the same 
behavior that was found in the scalar case |64| i.e. that the relativistic and non relativistic approaches do not coincide 
even when describing systems with zero binding energies as far as they interact with massive exchanges. 

The J = 1 state displays the same kind of departures from the scalar case than J = 0. Functions for a = 0, 1 
have been calculated using the values aps = 60, /i — 0.25 and A = 1.3. Contrary to the scalar case, binding energies 
are sizeably different: B a= o = 0.103 whereas B a= i = 0.0494. The physical wave function is obtained using the same 
procedure than for the scalar case, i.e. compute f^iQ and extract from them the coefficients Cj. Their values, 

Co = 0.749 and c\ = 0.662, are different from Co ~ c\ with Co bigger than c\. The averaged binding energy 

is B = 0.0793. The corresponding solutions are plotted in Figs. El One can see that /i dominates at small momenta 
(k << 1) but starting from k ~ 1, the components of relativistic origin become larger than f\. 

The splitting in binding energies is much bigger than for the scalar coupling. It can be seen in Figure ITBTi where 
the results of B a (a) for both A 2 eigenstates are plotted. The energy differences remain important even in the B — > 
limit - Figure [T5b - in accordance with the analytical considerations in section IVT1 

In summary, as it was noticed in Section IVII pseudoscalar coupling displays the largest deviations with respect to 
the non relativistic dynamics. Small and large spinor components are mixed to the first order. The coupling between 
fi and /2 is essential even for very weakly bound systems, the components of relativistic origin dominates already at 
moderates values of k and the splitting of the binding energies for the different projections of the J ^ states are of 
the same order than the energies themselves. 

IX. RESULTS FOR VECTOR COUPLING 

The stability analysis applied to vector kernels shows that vertex form factors are required for both J = + and 
J = 1 + states to obtain stable solutions. 

This is true in particular in the simplest application of vector coupling: the positronium J = 0~ state. The negative 
parity of the state comes from the intrinsic positron parity so that the corresponding kernels are those of the J 77 = + 
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two-fermion system already given in Appendix [X] In Table [I] are presented the values of the coupling constant a as 
a function of the sharp cut-off k max and for a fixed binding energy B — 0.0225. The dependence is very slow - 0.3% 
variation for k max S [10,300] - but it actually corresponds to a logarithmic divergence of a(k max ) as it can be seen 
in Fig. 1161 The origin of this instability is the coupling to the second component, whose kernel matrix element K22 
has an attractive, constant asymptotic limit. If one removes this component - which has a very small contribution in 
norm - calculations become stable and give for qnr — 0.30 the value ollfd — 0.3975. 



TABLE I: Coupling constant a as a function of the sharp cut-off k max for the J = positronium state with binding energy 
B = 0.0225 a.u. 
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The comparison of LFD ladder results with those obtained in perturbative QED or to the physical energies is 
meaningless due to the instability of the solutions themselves. The use of vertex form factors in a system of pointlikc 
particles would be hazardous and the introduction of renormalizable counterterms seems to be a more appropriate 
cure. 

First positronium results in Light Front Dynamics were obtained in |37l l38| . These authors introduced a large 
number of states in the Fock expansion but observed the same instability of the solutions. For a fixed value of the 
cut-off, the results become finite and can be compared. By taking k max = 10 and a = 0.3 - which corresponds to 
Bnr — 0.0225 - we found -Blfd = 0.0132, i.e. repulsive relativistic effects. The leading order QED corrections |79| 
reads 



Bqed - —r 



l + ^a 2 + o(a 4 ) 
16 



0.02516, 



and are so attractive. Equation (10) from [38| gives for k max — 10 the value Bjjlc — 0.0308, in qualitative agreement 
- thought still sizeably different - with Bqed- We should notice that a recent work [3^ analyzes the results of [3^ 
in terms of flow equations and obtains a closer value Bdlc = 0.02341. We conclude from that, that the ladder LFD 
predictions for such a genuine system are unable to reproduce even the sign of first order relativistic corrections. 
Because the lowest order corrections of the singlet state are not affected by the annihilation channels, the differences 
could be due to cross ladder graphes. 



For /1 ^ 0, the two fermion system is bound due to the /x-dependent terms (~ %Wy) in the vector kernel 1A8L 
since the /^-independent ones (Xij) are repulsive. This binding disappear in the non relativistic limit. 

When solving the equations for J = + state, the standard form factors l|18fl ~ depending on Q 2 and local in 
the non relativistic limit - were found to be insufficient for any power n to ensure stable solution. A Q 2 dependent 
gaussian form factor failed as well. This unstability comes from the /i-dependent terms. These are off-shell corrections 
depending on variables t, t' defined by 

Am 2 t = Ae 2 k - M 2 , 4mY = Ae 2 k , - M 2 (98) 

and are not regularized by a form factor depending on variable Q 2 . Such a function cuts off the high \k — k'\ 
components, but not the \k + k'\ ones. A similar situation is encountered in the framework of chiral perturbation 
theory [8(j and was solved by the replacement n(k, k') — > F(k)n(k, k')F(k'). 

Our way of doing is the following. Variable Q 2 — ~(k meson —ojti) 2 entering F(Q 2 ) is associated with the off-energy 
shell effects in the intermediate state containing one massive meson (p). In a similar way, we introduce the variable 
i] = m 2 — (ki — lot) 2 - see vertex 2 in the first graph of Fig. [31- and correspondingly 77' = m 2 — (k' 2 — lot') 2 from 
vertex 1. Variables 77, rj control the off-energy shell contribution to the fermion states and have been regularized by 
means of a cut-off function 



H(v) = 




This corresponds to a non-local form factor even in the non relativistic limit. On energy shell one has 77 = rj 1 = 0. 
Thus, for instance, the total form factor associated with vertex 2 in of Fig. |3| reads: 



Fm oc {Q 2 ,r 1 )=F{Q 2 )H{T 1 ) 



(99) 
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In center of mass variables Q the expressions for 77, 77' are: 



f(i_£«) 2m 2 i if - £s > 
I (1 + ^2) 2m 2 i if - £s < 



(100) 



and 



( 



(1 + ^) 2m 2 ?j' if 
(1 _ *L^») 2m 2 yj' if 



_ k£ > 




(101) 



with i,t' given defined in l|98l) . 

Each coupling constant is replaced by g — + gF(Q 2 )H(jj) - or 5 — * gF{Q 2 )H{rj') - and the kernel is multiplied 
by F 2 (Q 2 )H(r])H(rf). The values for A and n in H are taken the same than for F(Q 2 ), but could in principle be 
different. 

By means of <|99(1 . the solutions become stable but we notice that the use of only one kind of form factor is not 
enough to ensure the stability. Wave functions corresponding to /i = 0.15 obtained with n — 1 and A = 1.3 in 1991) 
are displayed in Fig. El Binding energy is B=0.0225 and a v = 1.485. They have normal behavior and one remarks 
sizeable relativistic component /2 starting from k 0.5 with a strong ^-dependence despite the small binding energy 
of the state. 

Let us now consider the J T = 1 + state. Solving the J T = l + ,a = equations with the F(Q 2 ) form factor only, 
leads to the same anomalies than for J v = + . With the non-local form factor the situations is regularised. With 
parameters B = 0.050, [i = 0.25, A = 1.3 and n = 1 for instance, one has a coupling constant a = 6.18 and a well 
behaved wave function. The same happens for the J w = 1 + , a = 1 state. When using, with the same parameters, the 
non-local form factor i|99[l . we get a = 6.01. 

The mass splitting between the two a = 0, 1 projections is shown in Fig. 1181 One first remark the striking behavior 
of a a (B) curves, i.e. larger binding energies correspond to smaller values of the coupling constant a. This fact - 
which takes place also for J — + states - is a consequence of the M 2 -dependence of the terms driving the vector 

kernel k^j in (|A8|) . Its contribution is large, because of /1 2 in denominator. Increasing the binding energy - i.e. 
decreasing M 2 - increases t, t' factors, and results into smaller values of a. When the M 2 dependence in t, t' kernel is 
frozen - setting e.g. M 2 — Am 2 - the usual a(B) variation is recovered (dotted curve in Figure fT^> . When including 
the full dynamics, both a a {B) curves get close each other in all the variation domain B = [0, 0.5], as it was the case 
in the scalar coupling. However due to their peculiar behavior - flat and almost parallel - the splitting in binding 
energies corresponding to a fixed value of the coupling constant, can be very large. One can also remark in Fig. 1181 
the different values of a a at B=0 despite the fact that the systems of equations for a=0 and a=l have - as in the 
scalar coupling - the same non- relativistic limit. This difference is due to the 1/ 'fj, 2 terms in the kernel. They are not 
relevant at the (l/m)° order but are crucial for binding a relativistic two-fermion system by vector exchange. For a 
fermion-antifermion system with massless exchange, e.g. positronium, the splitting at B=0 disappears. 



We have presented the explicitly covariant LFD solutions for the bound state of two fermion systems in the ladder 
approximation. A method for constructing non zero angular momentum states has been proposed and illustrated 
by numerical examples. It is based on satisfying the angular condition by a linear superposition of eigenstates of an 
operator commuting with the LFD ladder hamiltonian. 

We have separately examined the different types of OBE couplings and found very different behaviours concerning 
the stability of the solutions themselves and their relation with the corresponding non relativistic reductions. 

Scalar coupling (Yukawa model) is found to be stable without any kernel regularization for the = + state 
and coupling constants below some critical value a < a c — 3.72. For values above a c the system collapses. For 
J 11 = 1 + state the solutions of both a=0 and a=l projections are unstable. Their energy splitting is very small even 
for binding energies (B) of the same order than the constituents mass and vanishes at B=0. The physical solution, 
satisfying the angular condition, has been constructed by a suitable linear combination of a=0,l states. LFD binding 
energies are found to be close to those given by their non relativistic limit, even closer than the case of purely scalar 
particles (Wick Cutkosky model extended to ji 7^ 0) . The comparison with the non relativistic solutions shows always 
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repulsive effects. The LFD wavefunction is dominated by the components which has non-relativistic counterpart. 
Extra components of relativistic origin remain negligible even at large values of the relative momentum (k > m). 

Pseudoscalar coupling is also found stable for J 71 " = + state. It displays a very strong dependence of binding 
energies as a function of the coupling constant: they vary from B=0.001 to B=0.500 (in constituent mass units) while 
the coupling constant changes from a — 55.5 to a = 58.5. This dependence is due to the coupling to the wave function 
component of relativistic origin. Vertex form factors are required for J n = 1 + states. LFD solutions, obtained with 
regularized kernels, presents large deviations with respect to non relativistic case, even for weakly bound states, and 
display a big sensitivity to the cut-off parameters. The LFD wave function is dominated by relativistic component 
at relatively small momenta (k < m). The coupling between different components is strongly attractive and can 
compensate the repulsive effects observed in the Yukawa model. Thus, relativistic corrections can be attractive or 
repulsive depending on the quantum number of state, the value of the binding energy and even the mass fi of the 
exchanged meson. The energy splitting between different projections of J = 1 states is large and remains at B=0. 

Vector coupling presents the stronger anomalies. For [i = it has been applied to positronium + state. It is found 
to be unstable and, once regularized by means of sharp cut-off, the ladder approximation gives relativistic corrections 
of opposite sign compared to QED perturbative results. This failure shows the poorness of the ladder approximation 
in one of the rare cases in which it can be confronted to experimental results. For fj, ^ the LFD solutions collapse 
even using local cutoff form factors. The reason lies in the strong non-localities of the ^-dependent terms in the LFD 
kernel. These terms have their origin in the massive vector propagator and manifest as off-shell corrections of the 
/i — kernels. They have been regularized using appropriate vertex form factors. The J = 1 + state has thus been 
calculated. This state is not bound in the non relativistic limit and its existence in a relativistic approach is entirely 
due to the /x-dependent terms in the kernel The importance of this off-shell terms is thus dramatic. In particular their 
energy dependence generates a decreases of the binding energy as a function of the coupling constant, what questions 
the very meaning of the interaction strength. The a a (B) dependence for different projections of J = 1 states remain 
very close to each other even for B ~ to but their particular form - smooth and almost parallel variation - can give 
rise to large energy splitting for a fixed value of the coupling constant. 

Some general additional remarks concerning the relativistic calculations are given in order. 

(i) Contrary to the non relativistic case, vertex form factors are unavoidable in any realistic calculation. The full 
spinor structure generates highly singular kernels which are not regularized by local vertex form factors. It is 
clear that specially at large k-values, the obtained wave function and consequently the electromagnetic form 
factors will crucially depend on the way the regularization is performed. The large momentum components will 
thus be determined not by the dynamics but by uncontrolled parameters. We believe that here is the main 
drawback of relativistic approaches. 

(ii) The consequences of implementing the Lorentz invariance in a quantum mechanical description of a system 
are not only kinematical but mainly dynamical. Large differences with respect to the non relativistic solutions 
appear even in the zero binding limit for systems with ^ < 1 as far as the exchanged mass is non zero. We 
have explicitly shown for scalar and pseudoscalar couplings that the behavior of a(B) at B » differs from 
their non relativistic counterparts, a result already found in the Wick-Cutkosky model [64|. 

(Hi) The question about the sign of relativistic effects has no simple answer. They can be different, following: the 
nature of the constituents, the kind of interaction, the quantum numbers of the state, its binding energy, and 
even the mass of the exchanged particle. This shows that there are no simple recipes to perform a priori 
evaluations. 

(iv) The splitting of different projections of J=l states is very different following the kind of coupling. In nuclear 
physics - where the weight of scalar mesons in the binding energy is dominating - is expected to be very small. 
The same is true for the massless vector coupling like one-photon or one-gluon exchange. It can be however 
very large in relativistic models where pseudoscalar exchange plays an important role. 

Finally we would to emphasize one of the interest of using LFD in describing the relativistic composite systems. It 
lies in the fact that wave functions components appearing in this approach are closely related to their non relativistic 
counterparts. Some of these components are the formal equivalent of the usual non relativistic solutions while others 
are of pure relativistic origin. Relativity manifests both in modifying the former and in giving a sizeable weight to the 
latter ones. We have found that the coupling between these components plays an essential role, even in determining 
the stability of the solutions. In addition - except for the scalar exchange - the total wave function is dominated by 
the relativistic components at moderate values of its arguments (fc < to) and that, even for loosely bound systems. 
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APPENDIX A: KERNELS 



Kernels are obtained from equations (|45ll . 1461) . (|63|l and (|64|l as traces of 4x4 matrices. To calculate these traces, 
it is useful to express the scalar products between all the concerned four-vectors in terms of variables {k,k' ,6,8'). 
They read: 



kl 



to 
ki 

k[ 



k 2 = 

k[ = 

k 2 = 

k 2 = 



0, 

k~2 - — " k-^ - — - k- 

XLO-p, 

(1 - x)w-p, 

x' LU-p, 

(1 - x')uj-p, 

2el-m 2 , 



2e 2 k ,-m 2 , 



II 2 

— m , 



2e 2 (l - x) + -M 2 x, 



kvp 
k 2 -p 
k'rP 
k' 2 -p 

fci-fci = -kk 1 sin6>sin6>' cos0 + 2ef,,x + 2e k x' - 2e k xx' - 2s' 2 ,,: 



2e 2 x+-M 2 {l-x), 
2e 2 ,(l-x') + ht 2 x', 
2e 2 ,x + ±M 2 (l-x'), 



k 2 -k' 2 
k\ ' k<2 
k 2 -k[ 



-kk' sin6>sin6>' cos0 + 2e\x + 2£ 2 k ,x' - 2e\xx' - 2e\,xx' , 
kk' sin 9 sin^' cos0 + 2e 2 k (l - x)(l - x') + 2e 2 k ,xx' ', 
kk' sin 9 sin 9' cos <j> + 2e\, (1 - x) (1 - x) + 2e\xx' , 



where 



:(1 



kl-— cose') 

2 By 



(Al) 



(A2) 



Using the above result, we have obtained the analytical expressions of k,j kernels for J n = + , 1 + states. They are 
written below, coupling by coupling, in the form 



Kij(k,e,k',e',cp') = ClJ {k,e,k' ,6') + d. l0 {k,6,k' ,6')cosip' + e lj {k,6,k' ,9')cos 2 tp' 



(A3) 



with coefficients Cij,dij,eij invariant under the transformation (i, 
notations 



(i',k',6'). We introduce for shortness the 



s9 = sin 9 c9 — cos 9 S9 = k sin 9 C9 = k cos 9 cip' = cos ip' 

plus corresponding primed - and the following quantities 

m 2 (e 2 k +s 2 k ,)±2e 2 k e 2 k , 
e k ±m 

e k ^ 6 k' 



4 



A 



Coupling constants appear through a = 



1. Scalar 



Kernels for the scalar coupling were already given in 68] and are included here for completeness. 
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J = 0+: 



' ;: - -[m 2 A+ + 2e k e k >{e k e k , - C9C9')} + A+S9S9'ctp' 
-mA~ (S6' + S6c<p') 



an 

«12 



an 

— = +mA-(S9 + Se'ap') (A4) 
an 

K22 _ A+eaorf r™2, 

Q!7r 



A+S0S0' - [m 2 A+ + 2e fc £ fc ,(£ fc £ fc , - C9C9')]c<p' 



J = l+, a = 0: 



— = [2kk'e k e k > -b\c9c9'\ - s k e k ,A + s6s6'ap' (A5) 
an 

— = me k >(2el + A + )c6s6' - me k {2el, + A + )s9c9'cip' 
an 

— = me k (2e 2 k , + A + )s6c9' - ms k >(2s 2 k + A + )c6s6'ctp' 
an 



-£ k £ k '^ + s6s& + [2kk'e k e k > - b\c9c9'\ cip' 
J= l+,a= 1: 

— = - {me feS 2 0(A+ + 2e 2 ,) + m £fc , s 2 0'(A+ + 2e 2 ) + (c 2 + c 2 6')b 2 + - Ae k e k ,C6C6'} 
- {A + {elel,c6c6' - kk')s9s9'} off - { A + e k e k ,s 2 9s 2 9'} c 2 <p' 

CH7T 



= [me k (A+ + 2e 2 ,) - b 2 + ]s 2 9 - [me k ,(A+ + 2e\) - b\\s 2 9' 



- {kk'(e k - e k ,) 2 + (e k + £ k ,) 2 e k e k ,c9c9'} s9s9'c(p' + {e k ,{e k + e k >fc 2 9' - s+(e k - e k >) 2 } e k s 2 9c 2 <p' 

= -2e k e k ,C9S9' - ek r,{mA + -2e 2 k e k ,)c9's9' 

+ {2e k e y kCe' + e k c9 [e k ,{A+ - 2me k ) - e k ,{e k + Ewfc 2 9'\ } s9ctp' - (e k + e k/ ) 2 e k e k ,s 2 9s9'c9'c 2 tp' 

./o„. . 

: A" {(m + e^s 2 6»)S*6i' + mS'6'c^'-e A :s 2 6'5'6''c 2 ^'} 



an 

^ = [me k ,(A+ + 2e 2 )-b 2 + ] S 2 9' - [me k (A+ + 2e 2 ,) - b\]s 2 9 
an 

- {kk'(e k - eyf + (e k + e kl ) 2 e k 7,e k -c9c9'} s9s9 V + {s k (e k> + e k fc 2 9 - e+(e k > - e k ) 2 } e k ,s 2 9'c 2 ^ 
^ = m(e k + e k ,f - e k (mA+ - 2e k e 2 ,)c 2 9 - e k ,(mA+ - 2e 2 k e w )c 2 9' - Ae k e k ,C9C9' 

- (e k + e k ,) 2 (e k e k , - C9C9 r )s9s9 / c^ 

- {2b\{c 2 9 + c 2 9') + 2m[e fe (A+ + 2e 2 k ,)s 2 9 + e fc -(A+ + 2e 2 k )s 2 9'] + (e k + e k/ ) 2 e k e k ,s 2 9s 2 9 / - 8e k e k ,C9C9'} c 2 V ' 
{2k'e k e k ,C9 + e k ,(mA+ - 2e 2 k e k ,)c9'} s9' - {e k c9 [e k ,(e k + e k ,) 2 c 2 9' - e k ,(A+ - 2me k )\ - 2ke k e k ,C9'} s9cip' 



an 



- {e k ,c9' ((e k - e k ,fe+ - {e k + e k ,fe k c 2 9) + Ak'e k e k ,C9) s9'c 2 ip' 



V ^ 24 = +A-{(e k ~s k c 2 9)S9'- m S9 a p'-(s+-s k c 2 9)S9'c 2 ^} 



an 

= -2e k e k 'S0C0' -e7{mA+ - 2e k el)s9c9 
an fc 

+ s9' {s k c9' [(mA+ - 2e\e kl )c 2 9 + e k (A+ - 2ms k >)s 2 9] + 2e k e k .kC9'} op 1 - A + e k e k ,s9c9s 2 9'c 2 <p' 

= {2ke k e k/ C9 / + s k (mA+ - 2e k e 2 k ,)c9) s9 - {e k ,c9' [e k (e k + e k ,) 2 c 2 9 - e k (A+ - 2me k ,j\ - 2k'e k e k ,C9} s9'c V ' 

- {e k c9 [(e k - e k ,) 2 e+ - (e k + e k ,) 2 e k ,c 2 9'] + Ake k e k ,C9'} s9c 2 v ' 

— = {e k e k ,(A + -2me k >)c 2 9 / +e k 'e k (A + -2me k )c 2 9~e k e k ,(A + +2m 2 )- [(e k +e k >) 2 e k e k ,c9c9 , -2kk / e k e k >] c9c8'}ap' 

- [(e k + e k ,) 2 e k e kl c9c9 / - 2kk'e k e k >] s9s9'c 2 ip' 
^ = +A-e k s9c9S9'(l-c?<p') 
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^± - -A" {{m + e k , S 2 0')S9 + mS9'c V ' - e^ff SO<?ip'} 

= -A- {{e k , - e^,c 2 9 / )S9 - mS9'c<p' - (e+ - s^,c 2 9 / )S9c 2 ip / } 

^ = -A-ezse'ce'sen-c 2 ^) 

an fc 

— = [2e k e y {C6Ce' -e k e k >) - m 2 A+] ap' + A + S6S6'c 2 ip' 



J = 0+: 



Pseudoscalar 



— = -[m 2 A+ -2e k e k >(e k e k , -C6C6')} + A+S6S6'ap' 
an 

«12 



-mA~ (56*' - 56»V) 

+mA~ (S9 - S6'ap') (A6) 
A+SflSf?' + [m 2 A+ - 2e fc£fc '(£fc£fc' - C6C6')]ap' 



an 
W21 

Q!7r 
Q!7r 



J = 1, a = 



«12 

K21 



-{2kk'e k e k > + b 2 _c6c6') + e k e k >{A + - 2m 2 )s6s6' 'op' 
-mA" (e k >cQsQ' - e k s6c9'c(p') 

+mA~ (s k s6c6' - e k >c0stfap') (A7) 



an 

— = -e fe £fc'(A + - 2m 2 )s6s6' + {2kk'e k e k , + b 2 _c9c6')c(p' 
an 

J= l+,a= 1: 

H^il = 4C6»C6»'e fe e fc , + mA" (e fc s 2 - e fe / s 2 <9') + b 2 . (c 2 6> + c 2 6>') 
cot 

+ (efc - s k ,) 2 {kk' - e^s^cece')se S e'ap' + (e k - e k ,) 2 e- k E- k ,s 2 6s 2 6'c 2 tp' 
^ = (c 2 6 - c 2 6')b 2 _ - mA-(e k s 2 6 + e k> s 2 9') - {kk'(e k + e k ,) 2 + cOce'iek - e k ,) 2 e-e~,} sesO'c^ 

- {(b 2 _- me k A-)(l + c 2 6') + £ fc '£fe(A + + 2me k )s 2 6'} s 2 9c 2 (p' 

- {[MA - + 2ms k )-e k r,(e k - e k ,) 2 c 2 9'\e- k c9 - 2ke k e k ,C9'} s9c V ' 

- {(b 2 _+ me k ,A-)c6' + 2kk'e k s k ,c6} s9' + e k e k ,{e k - e k >f 's 2 9s9' ' c& 'cV 

\[2~K\i 



an 
2k 2 i 



A" {-{e k - e+c 2 9)S9' + mSOcip' + e k s 2 9S9'c 2 tp'} 



(c 2 6 - c 2 6')b 2 _ + mA-{e k s 2 e + e k ,s 2 9') - {kk'(e k + e k ,) 2 + c8c6'(e k - e k/ ) 2 e k s^,} s6s6' 'op' 

2 1 \-\(-\ 1 „2n\ 1 _ „-/A+ 1 1™,- \„2/)l J2nl„2,J 



- {(b 2 _+ me k ,A-){l + c 2 9) + e k e k ,{A+ + 2me k >)s 2 9} s 2 9'c 2 p' 

— = -mA-{e k s 2 6-e kl s 2 e')-b 2 _{c 2 6 + c 2 6')-Ae k e k ,C6C6' 
an 

+ {kk'(e k - M 2 + c9c9' [mA-(s k - e k >) -b 2 _- e k e k ,(A+ - 2m 2 )] } s0s6'ap' 

+ {%e k e k ,C6C6' + 2b 2 _{c 2 9 + c 2 9') + 2mA-{e k s 2 9 - e k ,s 2 6') + s 2 6s 2 6'(e k - e k/ ) 2 e k e w ) c 2 <p' 

^ K23 = \2k'e k e k ,C0+(b 2 _+me k/ A-)ce'}s8' 

an L ' 

+ { [M A+ + 2me k )s 2 e' + (mA+ + 2e k e 2 k ,)c 2 9'} s^c6 - 2ke k e k ,C9'} sOap' 

+ { [e k (A+ + 2me k ,)s 2 9 + (mA+ + 2e 2 k e k ,)(l + c 2 0)] e k ,c9' - Ak'e k e k >C9} sO'c 2 ip' 
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= -A" {{e k - c 2 9e^)S9' + mS9cip' - (e+ - e^c 2 d)S9'c 2 (p'} 

= { [e k (A+ + 2me k ,) - e^(e k - e kl ) 2 c 2 9] e^cff - 2k'e k e k ,C9) sO'ap' 
- {{b 2 _- me k A-)c6 + 2kk'e k e k ,c9'} s9 + e^e^ {e k - e k ,f s9c9s 2 9' 'c 2 tp' 

= {2ke k e k >C0' + (b 2 _-me k A-)c6} s0 

+ { [s k {A+ + 2ms k ,)s 2 9 + (mA+ + 2e 2 k e k <)c 2 9] s^,c9' - 2k'e k £ k <C9} s9'ctp' 

+ { [e k >(A + + 2me k )s 2 9' + (mA+ + 2e k e 2 k ,)(l + c 2 9')] 6^9 - Akk 1 'e k e k ,c9'} s9c 2 (p' 

= {-e k ,e k , [(A+ - 2m 2 )s 2 9s 2 9' + 2C9C9'} - b 2 _c 2 9c 2 9' + mA~ [e k ,s 2 9'c 2 9 - e k s 2 9c 2 9'\ } op' 

+ [2kk'e k e k , + e^e^,(e k - e k ,) 2 c9c9'\ s9s9'c 2 p' 

= -A-s-s9S9'c9s 2 <p' 

= -A- {-S9(e k , - e~c 2 9') + mS9'ap' + S9s 2 9' e'c 2 ^} 

= A- {{e k , - e k r,c 2 9')S9 + mS9'c<p' - (e+ - e- k ,c 2 9')S9c 2 ^} 
= A'e^s9'S9c9's 2 p' 

= (-b 2 _ - 2e k e k ,C9C9')ap' - A + S9S9 , c 2 ip / 

3. Pseudo- Vector 

Pseudovector kernels will be given as a sum of the pseudoscalar ones plus a term Sij which depends on variables 
t, t' defined in and vanishes on energy shell (t = t' = 0). 

K ij = K ij + $V 

The following expressions for Sij are valid only for x — x' > - with x,x' defined by (|A2|I - and because of that, 
coefficients <|A3ll are not symmetric in the exchange (i,k,9) <-> (i' ,k' ,9'). For x — x' < 0, the corresponding expressions 
are obtained by replacing t — > —t', if — > —t and their symmetry properties restaured. 

J = 0+: 

— = m 2 \m 2 tt' + (t-t')A- - (t + t')(e k >C6' - s k C9)\ - m 2 tt' S9S9'c<p' 
an 

c 

— = m \m 2 tt' + el(t ~ t 1 ) + e k {t + t')C6\ S9' + m \m 2 tt' - el(t - f ) - e k ,C9'{t + t 1 )} SOap 1 
an 

c 

— = m \m 2 tt' - el(t - t') - e k , (t + t')C9'\ S9 + m {m 2 tt' + e 2 Jt - t') + e k (t + t')C9\ S9'c<p' 
an 

c 

— = m 2 tt'S9S9' - m 2 \m 2 tt' + (t - t')A~ - (t + t')(e k ,C9' - e k C9)} ap' 
an 

4. Vector 

Vector kernels are written in the form 

o 

2 I m 

Kij = 2m tt — Vij + Xij 

in which Xij correspond to the /i = case. The contribution, due to /i-dependent term in the vector propagator, 

(PS) 

appears as being of shell corrections. Positronium kernels are simply given by — ~Xij- 



V2k 24: 

an 
an 

an 



^33 

an 

«34 
an 
-\/2k4i 
an 

an 

K43 

an 
K44 

an 
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J = 0+: 



£L = 2m 2 tt'^(m 2 + SeSe'ap') + (b 2 _ - 2e 2 ,e 2 ) 

2 

h ' L - - 2m 3 tt' 1 ^(S0' - S6cip')+mA-S6' 



2air /i 2 

= 2mht >™ 

2an /j, 

2a7r ii 



= 2m- i tt' — (S6-S6'c(p')~mA-S6 (A8) 



2m 2 tt' — (S9S9' + m 2 ap') - (A+S9S6' + 2e k e k , (e fe e fe < + C6C6')cip') 



J = l+,a = 0: 



2m 2 tt'^j [m 2 c9c9 / + e k e k , s9s9' cip'} + [m 2 A+c9c9 l + 2e k e k ,{kk' + 2m 2 sQsO' ap')} 



Kll 

2a-K ~ ' [i 

2 

= 2m 3 tt'^V Vk'cBsB' - e k s9c9'ap'] - me k , \A + c9s9' - 2e k e k , s9c9' VI 
2aw n z L J 

2 

= 2mht'^ T [e k s9c9' - e k ,c6s6'ap'} - me k \A+ s6c6' - 2e k e k ,c9s9'ap / ] (A9) 
2a7r /j, 

2 

-^1 = -2m 2 tt'^V |efcek'S0s0 / + m 2 c6»cf9'c^'l + \e k e k , A+s9s9' + 2t k t k ,{kk! + e k e k ,c9c9')cip'] 
2a-K [i A L J L J 

J = l+,a= 1 : 

The ^-independent kernels Xij are given by: 

^ = 2e fe£fe , [2C0C0' + (mA+ + e k ,e k c 2 9 + e k e w c^9 r )\ + [Akk'e k e k , + A+{kk' + e k el,)}s6s6' c^p 1 

+ 2e k e k >e k e k ,s 2 9s 2 9'c 2 tp' 

^ = [kk\4e k e k/ - A+) + e k e k , A+c0c0'] stfstf'cv' 

- 2e fe e fe .e A :(2e fe . - e£«V)« 2 0cV + 2s k e k , [e k £ k ,(c 2 9' - c 2 9) + m(e k s 2 9' - e k ,s 2 9)] 

= -2e k e k , {kk'cO + e k e k ,c9')s9' - {e k A+{m + e k ,s 2 9')c9 - Akk! e k e k ,c9'\ s6ap' + 2e k e k ,e k e k ,s 2 9s9 l 'cff 'c 2 <p' 

V2a7r 

Xu 
\f2aix 

X21 



= mA-S9c<p' 

-A+)- 

- 2e k e k <e k ,(2e k - e k s 2 9)s 2 9' c 2 tp' - 2e k e k > \_e k e k \(?9' - c 2 9) + m{e k s 2 9' - e k <s 2 9)] 



^ = [kk'{As k e kl - A+) + e k s k ,A+c9c9'] s9s9'ap' 



^ = -2e k e k , {m(e k + e k ,) + e k e k r,c 2 9' + e k ,e k c 2 9 + 2C9C9'} + {kk\Ae k e k > + A+) + e k e k ,A + c6ce'} s9s9' C ip' 



an 



+ 2e k e k , {(e+ + e k ,c 2 6'){e+ + e k c 2 9) + AC9C9'} <V 



U;; - 2e k e kl (kk'c9 + e y e k c9 r )s9' - 2s k e k , {2kk'c9 + (e+ + e- C 2 9) £ -,c9'} s9'c 2 ^' 



\f2aix 

+ {Akk'e k e k ,c0' - e^A+im + e w s 2 9')c9} s9ap' 
= n,A S'lc,:' 



\f2aix 

X31 

\/2a7r 

X32 



= -2e k E k > {kk'c9' + e k >e k c9)s9 - {e k ,A + (m + e k s 2 9)c,9' - Akk'e k s k >c9} s9'c^' + 2e k e k >e k e k ,c9s9s 2 0'c 2 (p / 

= 2e k e k ,{kk'c9' + e k e k ,c9)s9 - 2e k e k , {2kk'c9' + (e+ + e k ,c 2 9 / )e k c9} s9c 2 tp' 
V2a7r 

+ {Akk'e k e k ,c9 - e k ,A+{m + e k s 2 9)c9'} sO'ap 1 

^ = Ae k e k ,{kk' + e k -e k -,)s9s9'c 2 v l + 2{Ae k e k ,C9C9' + A+{e k e k , + e k e k ,c 2 9c 2 9 / - e k ,e k c 2 9 - e^c 2 ^)} 
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X34 = Q 

an 

= —mA~S9'ctp' 



— A - on' .J 



\p2~an 

/An 

= mA~S9cip 



^ = mA- -' 
an 

X43 = Q 

an 

^ = 2(244 -£W 

an 

and the Vij contribution reads: 

— = m[elc 2 e + e~ k ,c 2 9' - (e k + s k ')} - [kk' + e^s^c9c9']s9s9'ap' - £' k £~ k ,s 2 9s 2 9'c 2 <p' 



an 



V2< 



an 



an 



m[(e k - e*/) - e^c 2 + e^cV] + [kk' - e^e^c9c9']s9s9'c^ - e^(e+ - e^,c 2 9')s 2 9c 2 (p' 



-p^- = -nis^se'cff + E- k {m + El,s 2 e')s9c9cL P ' -els k r,s 2 9s9'c9'c 2 Lp' 
\/2an 

= (m + e k rs 2 9)S9' - mSBap' - s 2 9S9'c 2 ^ 

y2an 

= -m[(e* - e*/) + ej^cV - £- c 2 9] + [kk 1 - e^eyc9c9']s9s9'ccp' - e^,(e+ - e^c 2 9)s 2 9'c 2 ^ 

= m[(s k + e w ) - £l<?9 - £ k ,c 2 9'} - (kk' + eleyc9cB')s9s9'ap' - (2m + £~ k s 2 9)(2m + e^, s 2 & )c 2 ip' 

' - :! - me^sB'cB' + e^s9c9(m + e^, s 2 9' )ctp' - e^(e+ - e' k c 2 9)s9'c9'c 2 ip' 



\J~2an 

-^i- = ( m + e~s 2 9)Se' +mS9cip' - (2m + £ - s 2 9)S9'c 2 ip' 
\J2an 

= -me~ k s9c9 + e~, (m + e' k s 2 9)s9'c9'cip' - £ k e~,s 2 9' s9c9c 2 tp' 

\/2an 

1 iJ - me' k s9c9 + s^s9'c9'(m + £ k s 2 9)cip' - (4 - e w ^9')s9c9c 2 ^' 



V33 = -2(m + e k -,s 2 9')(m + e k rs 2 9)ap' -2e k -£ k r,c9c9's9s9'c 2 tp' 



' - 2e7:s9c9S9's 2 ip' 



an 

= (m + £^s 2 9')S9-mS9'cp' -e~ k ,S9s 2 9'c 2 <p' 

^ = (m + e k ~,s 2 9')S9 + mS9'cip' - (2m + e' k ,s 2 9')S9c 2 ip' 

^ = 2£-s9'c9'S9s 2 ^' 

an fe 



- -2(m 2 + S9S9'ap')c<p' 



l\7T 



APPENDIX B: RELATIONS BETWEEN THE COMPONENTS OF J = 1 STATE 



The wave function of the J — 1 state is represented in two forms: in the form (|49|) with the components tpi and in 
the form (|51J) with the components fi. The formulas expressing the components ipi in terms of the ft, in approximation 
M 2m, are given in Appendix C from Here we give these relations for arbitrary M. Note that ips and tps only 
differ relative to 1571, We denote below z — cos 9. 



<Pi 

<P2 



m 2 (2e k + m) m 2 R R 

^ h + ,e k (e k+m) ^ h ~ h + Zh - ^ Zk) 



^(V2f 1 -h-f 3 -2zU) , 
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g = V2(2e fc -M) 2 fc ^ (2e fc +m)(2 £fc -Af) 2 ^ | (4 £ 2 + 8Afe fc + M 2 )fc 
y " 16e|(£fc + m) 16e|fc "' " 16e|(efe +m) 

3M ^ _ z 2 (2 £fc -Af) 2 fc 2 \ h + V3M / + z 2 (2 £fc -M) 2 fc 2 \ ^ 



4fc V 12Me 2 k (e k + m) J J Ak \ 4Me 2 k (e k + m) 
3m \/3m 

1 /3~m 2 

^ 5 - 2\2kT k h ' 

(2e k -M)\ R (4 £ 2 +8M £fc + M 2 ) ^ 
^6 = — 5 V2/i - h + z/ 4 - V3z/ 6 * f 3 . (Bl) 

8m£fc 8m£fe 

The state with J = 1, a = 1 is determined by eq. (|67|l as a decomposition in four orthogonal spin structures . 
These four structures are expressed by eq. (16811 in terms of six structures Sjp, defined in l|48l) . with the coefficients 
hij given below. These coefficient are found as follows. We substitute the formulas l|66ll into I|B1|I . then eqs. (|B1(I - 

into l|49|) . In this way, the way function tfffi is expressed in terms of the four functions <?{ , i.e., obtains the form 
I[fi7jl. The coefficients at the front of gj 1 ' are the structures . Collecting these coefficients, we find in terms 
of six structures Sj^, in the form of eq. 168(1 with the following coefficients h^: 

_ Vim 2 _ Vim V3(£ fc - m)(4£ 2 + M 2 )z 

11 4£ fe (£ fe +m)' 12 4£ fc ' 13 16e 2 k k 

v / 3(4£ 2 + M 2 ) 



h\i — /li5 — 0, /li6 — 



8e k m 



, v / 3m 2 [£ fc (l-z 2 )+rn(l + z 2 )] v^m 

"21 — ; ^ 7 ' 1 22 — — - 



4£ fe fc 2 (l-z 2 ) ' zz 4£ fc 

. _ \/3(4£ 2 +Af 2 )[£ fe (l-z 2 )+m(l + z 2 )]z _ V3r, 

ft23 - i6ip(l - z 2 ) ' h2A ~ UY^z 2 )- 

^(4£ 2 +Af 2 )(l + z 2 ) 

<l25 — U, ^26 — 



8£fem(l — z 2 ) 

V3m 2 z _ n V3(e k -m)(4e 2 +M 2 )z 2 

A 3i — - — : : — 7= = i ft 32 — U, /I33 



2£ fe (£ fe +m)V2(l-z 2 ) " 8£ 2 fcV2(l-z 2 ) 

V3m _ V3(4£ 2 +A/ 2 )z 

"34 = — , , =, "35 — U, ft36 — 



fc^/2(l-z 2 )' ' 4£ fc m^/2(l-z 2 ) ! 



V3W 



2 



^41 = '142 = ^43 = ^44 = 0, /145 = ; — ====, h A& = 0. (B2) 

2£ fe fc v /2(l - z 2 ) 
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FIG. 9: g? = i t ... t 4 solutions for scalar coupling with a = 1.18, fi = 0.25 and sharp cutoff at k max — 10. Corresponding binding 
energy is B — 6.0498 




FIG. 11: Wave function components fi of the physical solutions (a) as a function of k at 6 = 30° and (b) their ^-dependence 
at fixed k = 1 value. Calculations are for the scalar coupling with a — 1.18, /i = 0.25 and sharp cut-off k ma x = 10. Binding 
energy is B=0.0501. 
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FIG. 12: Wave function components (in logarithmic scale) for J = + state with B=0.001, /j, = 0.15 obtained with pseudoscalar 
coupling and form factor A = 1.3. 




FIG. 13: B(a) for pseudoscalar coupling and J™ = + state with fi = 0.5 and two different form factors compared to non 
relativistic results 




FIG. 15: Splitting of the J=l solutions for pseudoscalar coupling. Results correspond to fi — 0.25 and A=1.3, n=l 
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FIG. 17: Wave functions fi for a J n — + state in the vector coupling with p = 0.15 and using the non-local form factor 19911 
with n=l and A = 1.3. The coupling constant is a = 1.485 and the binding energy B = 0.0225. 
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FIG. 18: Splitting of the J = 1 + solutions for vector coupling with /i = 0.25 and form factors l|99|l with A = 1.3 n=l. Dotted 
lines correspond to a fixed binding energy (B=0) in t, t' off-shell variables of kernel 1A8I 



